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1 Introduction 

The calculations of heavy quarkonium decay rates are among the earliest ap- 
plications of perturbative QCDtl and have been used to extract information on 
the QCD coupling at scales of the order of the heavy-quark mass. A consistent 
and rigorous framework for treating inclusive quarkonium decays has recently 
been developed, superseding the earlier non-relativistic potential model. The 
factorization approach is based on the use of non-relativistic QCD (NRQCD) 
to separate the short-distance physics of heavv-quark annihilation from the 
long-distance physics of bound-state dynamics □ The annihilation decay rate 
can be expressed as a sum of terms, each of which factors into a short-distance 
coefficient and a long-distance matrix element: 

T(H X) = Y,C{QQ[n] X)(H\O n \H). (1) 



The sum extends over all possible local 4-fermion operators O n in the effective 
Lagrangian of NRQCD that annihilate and create a heavy-quark-antiquark 
pair. The short-distance coefficients C are proportional to the rates for an 
on-shell QQ pair in a colour, spin and angular-momentum configuration n to 
annihilate into a given final state X, and can be calculated perturbatively in 
the strong coupling q; s (toq). The long-distance factors on the other hand are 
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well-defined matrix elements that can be evaluated numerically using lattice 
simulations of NRQCDEl For small average velocities v of the heavy quark in 
the quarkonium rest frame, each of the NRQCD matrix elements scales with a 
definite power of v and the general expression of Eq. (|l]) can be organized into 
an expansion in powers of the heavy-quark velocity. 

At leading order in the velocity expansion, the NRQCD description of S- 
wave quarkonium decays is equivalent to the non-relativistic potential model, 
where the non-perturbative dynamics is described by a single long-distance 
factor related to the bound state's wave function at the origin.| The NRQCD 
approach improves upon the potential model calculations for S'-wave quarko- 
nia by providing a rigorous non-perturbative definition of the long-distance 
factor, so this can be calculated using lattice simulations, and by allowing the 
systematic inclusion of relativistic corrections due to the motion of the heavy 
quarks in the bound state. 

Using the NRQCD factorization approach it is, in principle, possible to 
calculate the annihilation decay rates of heavy quarkonium from first princi- 
ples, the only inputs being the heavy-quark mass toq and the QCD coupling 
a s {mQ). Fixing the heavy-quark mass from, for example, sum rule calcula- 
tions for quarkonia, the analysis of quarkonium decay rates can provide useful 
information on the QCD coupling at the heavy-quark mass scale. 

For a reliable extraction of a s (rriQ), however, the theoretical uncertainties 
in the calculation of quarkonium decay rates have to be analysed carefully, and 
the effects of higher-order corrections in both the perturbative expansion and 
the velocity expansion have to be included. Radiative T decays are particularly 
well suited to study the impact of higher-order corrections, since not only the 
total decay rate but also the photon spectrum can be measured and compared 
with theory. 

2 The leading-order decay rate T — > 7 + X 

At leading order in the velocity expansion, radiative T decays proceed through 
the annihilation of a colour-sing let n = 3 S{ 1) bb pair. (We use spectroscopic 
notation with the superscript in brackets denoting the colour state.) The 
direct contribution, where the photon is radiated off a heavy quark, is shown 
in Fig. |l|(a). The corresponding leading-order decay rate is given by 

^In the case of P-wave quarkonia, the potential model calculations at O(ctf) encounter 
infrared divergences, which can not be factored into a single non-perturbative quantity. 
Within the NRQCD approach, this problem finds its natural solution since the infrared 
singularities are factored into a long-distance matrix element related to higher Fock-statc 
components of the quarkonium wave function. 
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(a) 



(b) 



Figure 1: Generic leading order Feynman diagrams for radiative T decay: (a) direct 
contribution, (b) fragmentation contribution. 
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where e& is the charge of the b quarkS Up to corrections of 0(v 4 ), the non- 
perturbative NRQCD matrix element is related to the T wave function at 
the origin through (T|d( 3 5i)|T) w 2N C \ ip(0) | 2 with N c = 3 the number 
of colours. More information about the decay dynamics is provided by the 
shape of the photon energy spectrum, which, at leading order in the velocity 
expansion and for the direct term, can be predicted within perturbation theory 
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where < x 7 — E^/m^ < li To very good accuracy, Eq. (Q) can be approxi- 
mated by a linear spectrum: dr^P/dar 7 2x 7 . 

It has been pointed out recently that there is an additional leading-order 
contribution to radiative T decays, where the photon is produced by-, frag- 
mentation, i.e. by collinear emission from light quarks, see Fig. ||(b)B The 
fragmentation contribution, although of 0(aa A s ) in the coupling constant, is 
enhanced by a double logarithmic singularity ~ In (m^/A 2 ) ~ 1/a 2 arising 
from the phase-space region where both the light-quark-antiquark splitting 
as well as the photon emission become collinear. The complete LO photon 
spectrum is thus given by 



dr LO (T -> jX) 



dX-y 



(4) 



where the first term on the right-hand side denotes the direct cross section, 
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Eq. (g), while the second term represents the contribution from gluon frag- 
mentation into a photonEl The fragmentation function is sensitive to non- 
perturbative effects and has to be extracted from experiment. The relative 
size of the direct component and the fragmentation contribution to the photon 
spectrum has been studied at leading order and is shown in Fig. Q The LO 
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Figure 2: The leading-order photon spectrum T — > *yXU Shown are the direct 
component (dashed curve), the fragmentation contribution (dot-dashed curve) and 



the sum (solid curve). Units (T|d( 5'i)|Y)/m| = 6/n have been used, and a s 



0.2. 



fragmentation contribution is important in the low-x 7 region, but suppressed 
with respect to the direct cross section for x 1 > 0.3.[] 

It is important to point out that the total decay width, including frag- 
mentation processes, T(T — » jX) = L dx 7 dr/d;r 7 , is not an infrared-safe 
observable. For x 7 — > 0, the fragmentation contribution rises like l/a; 7 , char- 
acteristic of the soft bremsstrahlung spectrum, and cannot be integrated down 
to x~ l = 0. Therefore, in perturbation theory, only the photon energy spectrum 
can be calculated for x 7 ^ 0, after collinear singularities have been factorized 
into the fragmentation functions. 

At present, only the region x~ l > 0.3 is accessible to accurate measure- 
ment, because of the strong background from tt — » 77 decays at low x 1 . A 
comparison between the theoretical prediction for the shape of the photon en- 
ergy spectrum and the available experimental data can thus be restricted to 
the direct term. Such a comparison is shown in Fig. [|3 The theoretical curve 
has been modified to account for the experimental efficiency and energy res- 
olution and has been normalized to the data. The steep rise towards large 

c At next-to-leading order, the fragmentation contribution involves also quark fragmentation, 
which is much harder than gluon fragmentation and may be important at larger values of 
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Figure 3: Experimental data for the photon spectrum compared to the emcieriisy- 
attenuated and energy-smeared prediction for the direct process at leading orderu 



x 7 predicted by LO perturbative QCD is not supported by the experimen- 
tal data. In particular for x 1 > 0.7 the discrepancy is significant. Near the 
phase-space boundary i 7 -t 1, hadronization effects may become important, 
but they are not expected to modify the energy distribution significantly in 
the intermediate region a; 7 ~ 0.7. To accommodate the data, different models 
have been proposed, including a parton shower Monte Carlo approach, also 
shown in Fig. ^, or parpametrizations of non-perturbative effects in terms of 
an effective gluon mass™ A proper theoretical understanding of the photon 
spectrum, however, requires a consistent inclusion of higher-order terms in the 
perturbative expansion as well as a careful examination of the contributions 
from relativistic corrections. 

Given the above-mentioned theoretical problems, the present extractions 
of a s (rrib) from the total radiative decay rate have to be considered uncertain. 
Even neglecting fragmentation contributions at low cc 7 , model approaches have 
to be used to extrapolate the experimental data into the lower energy region to 
obtain the total decay rate, introducing uncontrolled theoretical uncertainties. 
Analysing the effect of higher-order corrections to the photon energy spectrum 
is thus not only an interesting problem of perturbative QCD per se, but will 
also lead to a more reliable determination of the strong coupling constant from 
radiative T decays. 
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3 Higher-order corrections to T — > 7 + X 

The general factorization formula for the T decay rate of Eq. (|l|) is a double 
power-series expansion in the strong coupling a s and the heavy-quark velocity 
v. Higher-order contributions to the LO estimate include both 0(a s ) correc- 
tions to the short-distance cross section and relativistic corrections 0(v 2 ) due 
to the motion of the b quarks in the T bound state. Subsequently, we will dis- 
cuss the different sources of higher-order corrections and examine their impact 
on the photon energy spectrum. 

3.1 Soft-gluon resummation 

It has been argued that potentially large logarithms ln(l — x 7 ) associated with 
the imperfect cancellation between real and virtual emission_of soft gluons as 
x 1 — > 1 may contribute to all orders in perturbation theoryEj The resumma- 
tion of these soft-gluon effects may then give rise to a Sudakov suppression 
~ 1/ exp(a s ln 2 (l — x~)) near the endpoint of the photon energy spectrum. 
However, a recent analysis reveals that the logarithms of (1 — x~) cancel at 
each order in the perturbative expansion, in the non-relativisticJimit where T 
decays proceed via the annihilation of a colour-singlet bb pairtU Accordingly, 
no Sudakov suppression arises at leading order in the velocity expansion. It 
has also been shown that the parton shower Monte Carlo approachp used in 
the recent determination of a s from radiative X decays, does not correctly take 
into account gluon radiation in higher ordersuil 

3.2 Relativistic corrections 

Although the velocity of the b quarks in the T bound state is small, v 2 rj 
0.1, relativistic corrections may contribute significantly in certain regions of 
phase space. The NRQCD approach allows a systematic calculation of these 
corrections, introducing however additional non-perturbative matrix elements. 
Thus, in the absence of lattice results, one has to resort to phenomenological 
inputs to estimate the size of the higher-order terms in the velocity expansion. 

The leading relativistic correction to the colour-singlet decay rate is of 
0{v 2 ) and pan be expressed in terms of the 'binding energy' of the quarkonium 
Mx — 2m;, La Depending on the value of the 6-quark pole mass, the 0{v 2 ) cor- 
rection may be as large as 2ja2o, but it affects the shape of the photon spectrum 
only in the region x 7 > 0.8113 At 0(v 4 ) in the velocity expansion, radiative T 

decays can proceed through the annihilation of colour-octet n — 1 Sq 8 \ 3 Pq 8 2 
bb pairs into 73 final states at 0(aa s ). The enhancement of the colour-octet 
short-distance annihilation cross section ~ nja s partly compensates the strong 
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suppression of the long-distance matrix element ~ u 4 . Since the LO short- 
distance process for colour-octet terms is kinematically restricted ~ 6(1 — x 7 ), 
colour-octet contributions to the direct component of the decay spectrum are 
important only in the endpoint region x 7 > 1 — v 2 ~ 0.9, taking into account 
an effective smearing of the delta function due to the resummation of higher- 
order terms in the velocity expansionEiLColour-octet processes also contribute 
to radiative T decays by fragmentationlij Both gluon and quark fragmentation 
turn out to be sizeable at low photon energies x 7 < 0.3. A separate measure- 
ment of the gluon fragmentation functions from the lower part of the energy 
spectrum, as suggested by the analysis of the colour-singlet decay channel, may 
thus not be feasible. 

To summarize, relativistic corrections to the photon energy spectrum, al- 
though potentially sizeable, are concentrated in the upper and lower endpoint 
region and do not significantly modify the shape of the spectrum in the inter- 
mediate energy range 0.4 < x 7 < 0.9. 

3.3 Next-to-leading order QCD corrections 

For intermediate photon energies, the spectrum is expected to be well described 
by the direct colour-singlet channel. The next-to-leading order corrections to 
the total decay rate of Eq. (||) are potentially large, depending on the choice 
of the renormalization scale. An earlier calculation yields 

rS°(T^7*)=r£ (i_ 21^(1.67 ±0.36)), M = 2m b , (5) 

with a large theoretical uncertainty of ~ ±20% coming from the numerical 
evaluation of the loop integralsll3 

The calculation of the higher-order QCELcorrections to the photon energy 
spectrum has been completed only recentlyO Generic diagrams that build up 
the decay rate in next-to-leading order are depicted in Fig. [|. Besides the usual 
self-energy and vertex corrections for photon and gluons (a) , one encounters a 
large number of box diagrams (b) as well as gluon radiation off heavy quarks 
(c) and gluon splitting into gluons and light-quark-antiquark pairs (d). The 
evaluation of these amplitudes has been performed in the Feynman gauge and 
dimensional regularization has been used to calculate the singular parts of 
the amplitude. The masses of the nif = 4 light quarks have been neglected 
while the mass parameter of the b quark has been defined on-shell. The ex- 
change of Coulombic gluons in diagram (4b) leads to a singularity ~ 7t 2 /2/3r, 
which can be isolated by introducing a small relative quark velocity (3r. The 
Coulomb-singular part can be associated with the interquark potential of the 
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Figure 4: Generic next-to-leading order Feynman diagrams for direct radiative T 
decay. 

bound state and has to be factored into the non-perturbative NRQCD matrix 
element. Only the exchange of transverse gluons contributes to the next-to- 
leading order expression for the short-distance annihilation rate. The infrared 
singularities cancel when the emission of soft and collinear_final-state gluons 
and light quarks, described by universal splitting functionsEJ, is added to the 
virtual corrections. The analytical result for the matrix element squared has 
been implemented in a Monte Carlo integration program, so that not only the 
inclusive decay rate but any distribution can be generated. 
For the total width we obtain at next-to-leading order: 

rg£°(T-> 7 jr) = iS£ [1 - ^ (i.n + AMin^)l , (6) 

where f3o{n\i) = (lliVc - 2nif )/30 The theoretical uncertainty due to the 
numerical phase-space integration has been reduced to < 0.5%. The new, moxe 
accurate result for the total width is consistent with the previous calculationEj 
within the estimated numerical error. 

The next-to-leading order corrections significantly flatten and deplete the 
photon energy spectrum for z 7 > 0.75 as shown in Fig. [5], indicating that 
the discrepancy between theory and experiment at large photon energies can 
be reduced by the inclusion of higher-order QCD corrections. A more detailed 
discussion of the photon spectrum at NLO and a comparison with experimental 
data will be presented in a forthcoming publication. 
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Figure 5: Photon energy spectrum for direct radiative T decays in leading and 
next-to-leading order (a s = 0.2). 

4 Conclusions and outlook 

The photon spectrum in radiative T decays is a very interesting laboratory to 
study perturbative and non-perturbative QCD effects. In the low-x 7 region, 
fragmentation contributions dominate, while near the upper endpoint of the 
spectrum rclativistic corrections and hadronization effects become important. 
The next-to-leading order results presented here indicate that the intermediate 
region of the energy spectrum can be accounted for by perturbative QCD, once 
higher-order corrections are taken into account. This should allow a more 
reliable extraction of a s from radiative T decays by restricting the analysis to 
the region of the energy spectrum that can be described by NLO perturbation 
theory. 
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